Unbalanced Renormalization of Tunneling 
in MOSFET-type Structures in Strong High-Frequency Electric Fields 
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Two-dimensional electron gas coupled to adjacent impurity sites in high-frequency out-of-plane 
ac control electric field is investigated. Modification of tunneling rates as a function of the field am- 
plitude is calculated. Nonlinear dependence on the ac field strength is reported for the conductivity 
of two-dimensional electron gas. It develops a periodic peak structure. 
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I. INTRODUCTION 

Rapid development of experimental techniques at 
nanoscalei^i^i^ has stimulated theoretical advances in 
describing quantum phenomena for various geometries 
and settings. Extensive study has been done on nonlinear 
effects in a few state quantum system subject to strong 
harmonic control^M^ . 11 . 12 ! 13 ^^ 15 ! 16 ! 17 . 18 ^ 9 such as a 
double quantum dotji 7 - an array of coupled quantum 
dotsji 8 - super laticeSfi^ etc. In this paper, we investigate 
the influence of the ac field on one-state quantum objects 
coupled to two-dimensional electron gas (2DEG) via tun- 
neling. 

The systems of such geometry have been recently 
used in experimental as well as theoretical study of 
few and single electron spin manipulations^ spin- 
to-charge conversion measurements, 4 Rudcrman-Kittel- 
Kasuya-Yosida (RKKY), 3 and quantum Hall 20 effects. In 
the mentioned experiments, 2DEG is formed by confine- 
ment in a two-dimensional layer grown in a pile between 
between the layers of the wide gap material. The elec- 
tron concentration in 2DEG can be varied significantly 
and is usually controlled electrostatically by split-gate 
technique. A similar system can be also created in the 
inversion layer in metal oxide semiconductor field effect 
transistors (MOSFETs). In both cases, the impurity cen- 
ters (sites), localized usually outside the 2DEG in adja- 
cent layers, play an important rolej 21 > 22 One of the first 
experimental evidences here was the observation of ran- 
dom telegraph noise in conduction of the inversion layer 
in MOSFET^i 

Surprisingly, little attention has been given to the con- 
trol of the impurity states, and thus the properties of 
2DEG, dynamically. Unlike the well-known phenomena 
of dynamical control of tunnelin g 13 ' 14 i 15 ' 16 i 17 ' 18 in few 
state electron systems, e.g., double quantum dot, the 
impurity-2DEG system provides more degrees of freedom 
to change properties and correlations which are not re- 
lated to tunneling directly. One of the examples here is 
the possible indirect influence over the RKKY interaction 
mediated by 2DEG electrons^ 
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In what follows, we demonstrate that periodic high fre- 
quency potential (electric field) applied perpendicular to 
2DEG leads to nontrivial renormalization and disbalance 
of the tunneling between the impurity sites and 2DEG. 
Moreover, tunneling modification, as well as Coulomb 
activation of the impurity sites, induces oscillatory be- 
havior of 2DEG conductivity as a function of the ampli- 
tude of applied periodic field. This variation is similar 
to Shubnikov-de Haas oscillations^ 3 - but have different 
underlying physics. 

In the next section, we formulate the impurity-2DEG 
model. Section IIIII is devoted to the construction of the 
corresponding stationary many-body problem using Flo- 
quet states. Time-averaged quantities of interest are de- 
fined. In Sees. IIVI and [Vl nonlinear dependence of tun- 
neling rates is obtained, starting with the simpler case 
that neglects 2DEG electron scattering on impurity. The 
scattering dynamics is analyzed. Finally, in Sec. lVIl field 
amplitude dependence of conductivity is found. This de- 
pendence, together with the expression of the tunneling 
rates, is the main result of the paper. 



II. MODEL 

As mentioned above, we consider zero temperature 
2DEG interacting with an impurity electron localized in 
the adjacent layer of a wide gap material. Both systems 
are subject to external plane polarized harmonic field, 
E z (t), with the frequency u>o and polarization along z 
axis, perpendicular to the 2DEG plane, see Fig. [T] The 
field is treated in the dipole approximation appropriate 
as soon as the characteristic size of the nanostructure is 
smaller compared to the field wavelength. The impurity- 
2DEG interaction is via tunneling between 2DEG and 
impurity states, as well as through Coulomb scattering of 
conduction electrons on empty, positively charged, donor 
impurity site. The electrons are considered spinless, and 
only the ground state of the impurity site is taken into ac- 
count. In many cases, the heterostructure has more than 
one impurity site next to the conduction layer. This sit- 
uation, including the effect of spins, is discussed in the 
last sections, and in most cases, the effect of the exter- 
nal field can be deduced from the one-site spinless model 
investigated below. 
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The unperturbed Hamiltonian is 
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(2.1) 



where Ud(r) is the impurity localization potential, and 
U2deg( z ) denotes the potential profile that forms 2DEG. 
The 8 functions are defined so that Hq i resemble the 
unperturbed Hamiltonian of an electron sitting on the 
impurity (8 i>c t = 1, ^deg = 0) and in the 2DEG (8^ d = 
0, <$i.2DEG — !)• The interaction between electrons as 
well as due to the external field is 
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where E z (t) = E z coswrj^ and q s represents the screening 
wave vector. Following the standard procedure, we define 
the amplitudes 
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and 



V kk > = -2 j drdr'|Vo(r)| 
x ^(r')^'(r') 



^fc(r) (2.3) 
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Here, ipo( r ) is the wave function of the noninteracting 
ground state localized on the impurity, i.e., ■ffo,dV'o( r ) = 
E i(j (r), and ^fc(r) corresponds to the fc-th state in 
2DEG, i.e., i?o,2DEG^fc(r) = E k ij: k (v). The shape of 
ipo( r ) and ^fc(r) depends on the form of U d {v) and 
C yr 2DEc( r )- Though it is not, strictly speaking, necessary, 
we will assume Hok to be independent of k and real, i.e., 
H k0 = Hok — > A, to simplify notations. The electron- 
electron interaction in 2DEG is ignored. With the above 
notations, we arrive to the Hamiltonian 

H = [E + Voo(t)}dU + J2[Ek + V 2 DEG(t)}c{c k 

k 

+ d^J2 V kk'clck>+Aj2(^ c k+cldy (2.5) 

kk' k 

where Voo(t) = —eE z (t) J drtpQ(r)zip ( r ) and, similarly, 
V2DEg(*) ~ / dTip k (r)zt/jf.( r )' The latter is assumed to 
be independent of k since only the distribution along 
z of the lowest band is of interest. The z distribution 
of |V>o( r )| 2 and |V>fc(r)| 2 is concentrated around, respec- 
tively, the impurity center and the middle of the L^deg 
quantum well. Therefore, zo,o = / drV , o( r ) z V , o(r) and 
z k,k = J dripKrfzipki'r) refer to the position of the impu- 
rity and 2DEG along the z axis. As we will see later, only 
the difference between these two quantities is of interest. 

One can adopt the following values for order of mag- 
nitude estimates. With the separation between the im- 
purity and 2DEG of the order of several angstroms, and 
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FIG. 1: Impurity site coupled to the 2DEG conduction elec- 
trons. Strong ac field of frequency coo is applied perpendicular 
to 2DEG. The ladder of quasienergies develops. All quasiener- 
gies above Ef and within the range of ~ Xujq are used for 
tunneling out process. The quasienergies below the bottom 
of the conduction band, E c , are not active. 



the barrier height of the order of several eV, the tunnel- 
ing amplitude will vary by ~ 1 10 meV and smaller 
depending on the distance. The frequency luq is of the 
order ~ 10 meV; the temperatures T < 1 K; the size 
quantization > 100 meV. These values are feasible for 
Si/Si02 structures. 



III. FLOQUET STATES 



The Hamiltonian (|2.5[) is periodic in time with the pe- 
riod 2tt/ojq. It is natural to utilize this symmetr y 24 ' 25 i 26 
Similar to space-periodic solid state lattice structures, it 
was showr* 2 ^ that the wave function corresponding to the 
periodic Hamiltonian is of the form 



\u e (t)), 



(3d) 



where \u £ (t)} — \u £ (t + 2tt / ujq)) , and e is the 
quasienergy. Moreover, it was demonstrated that the 
set of quasienergy states can be treated similarly to the 
conventional system of stationary eigenstates — i.e., the 
system initially set up in a certain quasienergy state (or 
distribution over these states) remains at the same state 
(or with the same distribution) over the entire evolution 
of the system^ 

The transitions between the quasienergy states corre- 
spond to the perturbations which break the periodicity. 
This has been investigated in the literature^ In our case, 
we assume the time scale of such perturbations due to 
environment (or other factors) to be much larger than 
the one of interest. As a result, one can analyze the 
quasienergy spectrum of the model to obtain the infor- 
mation about tunneling effects in the system. 

Let us show a few steps to support the above 
statement. It is convenient to use the interac- 
tion representation, factoring out the evolution due 
to the oscillatory part of the Hamiltonian (|2.5p . 
The corresponding evolution operator, Uo(t) = 

exp{ T ^sinuj t[Vood ii d + V2DEGi r 2 k CkCk}}, is still peri- 
odic so that one can define 

\Mt)) = u f (t)\Mt)) 



e iet \u s (t)), (3.2) 
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with \u £ (t)) = \u £ (t + 2tt/(jJq)). The corresponding 
Schrodinger equation is 

ij t \Mt))=H(t)\Mt)), (3.3) 

where 

H(t) = E dU + J2 E kcl c k+dSJ2 V kk'clck' 

k kk' 

+ (^c k e-* xanUot +4de 4Asin " ot ).(3.4) 

fe 

The effective strength of the external periodic field is de- 
fined as A = (Vbo — V2deg)/wo with Vbo and V2DEG rep- 
resenting the amplitudes of Vbo (t) and V2DEG [t) respec- 
tively. Note that Vbo and V2DEG have opposite signs and, 
thus, A is proportional to the average distance |^o.o — z k.k\- 
Our goal is to obtain the equation for quasicnergy. One 
can easily form equations for the Fourier transform of the 
periodic part of the quasicnergy wave function, \u e (t)} ~ 
J2 m e imuot \u e m ). From Eq. ^) we have 

J2 e- lm " ot [e + com H(t)} \uLJ = 0. (3.5) 

m 

Let us now recall that the time-dependent exponentials 
in Hamiltonian (|3.4p have a simple series representation 
el Asin^t = J]^ = _ oo J n (A)e mwot in terms of the Bessel 
functions. It should be noted that the harmonic form 
of the external field, and thus the time-dependent expo- 
nential in Eq. (|3.4|) . is not necessary. For any periodic 
zero-average field, one can define the above series repre- 
sentation. In this case, the Bessel functions J„(A) are re- 
placed with the coefficients, J n (A) — * / n (A), which carry 
the structure of a single oscillation. The results obtained 
below will be qualitatively the same with this modifi- 
cation. Since the harmonic field provides more insight 
into the physics of the phenomenon, we use it in further 
derivations instead of a more complicated time-periodic 
potential. 

Defining a vector column of the states as \v) = 
{..., |-«o), ...} T , we finally obtain the time- 

independent Schrodinger equation for quasienergies, 

H\v)=e\v), (3.6) 

where the stationary Hamiltonian is 

H = Epd^d + E k c k c k + loqIz + dS Vkk<c\c k > 

k kk' 

+ y, AJ «( A ) { l( - )c U + 4 n)rftc fc) • ( 3 - 7 ) 

Here the additional operators are understood, if one de- 
fines a column vector, |e m ), with all entries zeros except 
for the m-th entry which is "1." Then, I z \e m ) — m\e m ), 
I±\e m ) = |e m ±i). We use the superscript in parenthe- 
ses to generalize the power as Zj_ "= 1^ ' = I™ . The 



quasienergy spectrum is now a solution to the Kondo- 
type spin-assisted tunneling problem, where operators I± 
correspond to renormalized rising (lowering) operators of 
a large integer spin (S — > 00), or an asymptotically large 
ensemble of identical two-state systems. Note that I±, 
as they have been introduced to rewrite Eq. (|3.5p in the 
form (|3.6|) , are not quite the spin rising (lowering) opera- 
tors. Nevertheless, in the limit of large spin (5 — » 00) and 
finite magnetization, they differ only by a constant factor 
which has no effect on the subsequent calculations^ 

We should also demonstrate that the stationary prob- 
lem with Hamiltonian (|3.7p is sufficient to compute phys- 
ically observable quantities of interest. In this paper, 
we are after the tunneling process in between the impu- 
rity and 2DEG, as well as the conductivity of 2DEG, 
therefore, it is natural to investigate the dynamics of 
the average occupation number for the impurity site, i.e., 

(tp(t)\ddi \ip(t)) , or the amplitude of 2DEG electron tran- 

1 1 

sitions for states k, k , i.e., (ip(t)\c k c\,\il>(t)) . The time- 
average is over the period 2tt/lu of the fast external field 
oscillations. 

Taking into account the properties of the quasienergy 
states, we can focus on the average over a single 
quasienergy state. As mentioned above, the average of 
two conjugate operators of the same type is sought. This 
simplifies the expression further as (ip e (t)\dd^\ip e (t)) = 
(ipe(t)\dd^ \ip e (t)} . The time-averaged quantity becomes 

<^ e (t)|ddt|^ e (*))' = 5^<ti TO |5(0,t)d(t)^(t)5(t,0)|ti I n) 

rn 

(3.8) 

Here, the operators are in the interaction picture and 
evolve according to the first three (main) terms of the 
Hamiltonian (|3.7p . while the standard scattering ma- 
trix is due to the perturbation — the last two terms in 
Eq. (|3Jl) . 

As a result, the dynamics of the average occupation 
probability at the impurity site is entirely determined 
by time-independent Hamiltonian (|3.7p . Similar argu- 
ments hold for the transition amplitudes in 2DEG and 
thus the conductivity. A standard equilibrium proce- 
dure of switching the interaction "on" adiabatically from 
t = —00 can be used. In this case the initial dy- 
namics is stationary in the first place, \ip e (t)) = const, 
since Hamiltonian (|3.4p becomes time- independent. This 
makes \u m ) = \uo)^o,m- The expression for the average 
becomes 

MS))' = (3-9) 
= (Go, e |5(-oo, t)d(t)d'(t)S(t, -oo)|w , e ), 

where \uq) is the usual initial state for nonintcr- 
acting fermions. In what follows, we will use 
the shorthand notation for the complete average 
(S(-oo,t)d(t)d' ! '(t)S(t, -00)) instead of the one m 
Eq. jSU). 
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IV. NONLINEAR TUNNELING WITHOUT 
SCATTERING 

In this section, we obtain the tunneling rates consid- 
ering Hamiltonian (|3 . 7|) without Coulomb scattering on 
the impurity, i.e., the fourth term. Let us explicitly show 
the main part, 

H = E dU + E kc k c k + oJoh (4.1) 

k 

and the perturbation, 

V = A JnW (l ( - ] c\d + I [ + ] Sc k ) . (4.2) 

n.k 

The perturbation (|4.2[) leads to equilibration of the im- 
purity occupation probability P. Using this state for av- 
eraging, we have (S(— oo, oo)) = e 1 ^ . As a result, the 
tunneling rates can be found by calculating the zero- 
temperature impurity electron self-energy. The Green's 
function G{t,t') = -i{Td(t)d) (t')S(oc, -oo)) is of inter- 
est, where Sfah) = Texp[-i J t * 2 dt'V(t')}. 

For small A, one-loop approximation is sufficient^ 
In Appendix [A] we calculate the self-energy for higher 
orders in A. However, they do not introduce any new 
physics and may be omitted. The tunneling rate 7 is 
given by imaginary part of the self-energy, 

£iH = £A 2 .£(A)<7fc(w-mwo), (4.3) 

where gk(u) is noninteracting Green's function of 2DEG 
electrons. The tunneling in and out from the impurity 
can be clearly separated. The result is 

M? n/out 

lin/out = -2ttA 2 D 2 J ™W- ( 4 ' 4 ) 

Here, D 2 is the density of states for 2DEG. The limits are 
M} n = -0(Eo-E F ), Ml = - YTn=i F 0(E -nuj -E F ), 
and Ml ut = 6{E F - E ) + Mf n , M 2 ut = 00, where 
Ep is the Fermi energy of 2DEG. These limits and the 
summation terms have a clear physical meaning. They 
correspond to tunneling in and out from the impurity 
quasienergy states, see Fig. [TJ At the same time, they 
may be viewed in terms of allowed multiphoton processes 
from the 2DEG state below the Fermi surface (tunnel- 
ing in, negative m) and to the empty states above Ep 
(tunneling out, positive m). However, the latter lan- 
guage should be used keeping in mind the explanation 
via the quasienergies. The actual photons also account 
for the renormalization of the tunneling amplitude, which 
is done automatically in our treatment. The infinity in 
M 2 ut is true in the limit sense, lim^^^iSoo/ Auio, and 
denotes the upper edge, Eoo, of the 2DEG band or the 
next conduction band if present relative to the external 
field strength Xuq. 



To be specific, let us discuss the case when the electron 
concentration in the conduction channel, no, is low, such 
that the Fermi energy of 2DEG (with respect to the bot- 
tom of the conduction band, E c ) is smaller as compared 
to the external field frequency. The magnitude of the 
field exceeds the latter and is much smaller than E^. 
The chain inequality is Eq ~ E F < u>q < u>o A <C £oc, 
where all the energies are measured from the bottom of 
the 2DEG conduction band, E c . This is a natural as- 
sumption for many 2DEG systems used for few electron 
manipulation in recent experiments. The tunneling rates 
are 

7i„ = -27rA 2 £> 2 J 2 (\)6(E F - E ) (4.5) 

and 

00 

lout = -2nA 2 D 2 j2 mW- ( 4 - 6 ) 

m=6(E F -E a ) 

For a weak external harmonic field, A <C 1, both 
tunneling rates approach the well known result 70 = 
—2tt A 2 D 2 9(±E F =F E ). This is also true for the more 
general case given in Eq. (|4.4[) . The low amplitude of the 
field suppresses the multiphoton absorption (emission) 
process, J^yiW ~ * Oj and allows the single photon pro- 
cesses as a small first-order correction, J 2 (A), while the 
renormalization of the non-assisted tunneling vanishes, 
J 2 (A) - 1. 

Larger external field amplitudes activate more 
quasienergy states. Keeping in mind the earlier discus- 
sion, this process can also be viewed in terms of induc- 
tion of multiphoton transitions with the maximal num- 
ber of photons absorbed (emitted) per transition ~ A. 
In this case, the tunneling rates depart from each other, 
see Fig. [51 In the limiting case of A ^> 1, the tunneling 
out rate becomes 7 ou t — > —2itA 2 D 2 \ 1 while ji n oscil- 
lates according to Jq(A) and converges to zero as ~ 1/A. 
Similar results can be obtained directly by averaging the 
transition matrix element due to Hamiltonian (|2.5[) over 
the fast oscillations of driving field, keeping the terms 
of the order A 2 , see Appendix [Bl The above approach, 
however, provides more physical insight and is more con- 
venient for further discussion. 

It should be noted that the above results apply in a 
more general case when Eq is significantly higher or lower 
then E F . In this case one can define Eq = Eq + fhujQ, 
such that \Eq — Ep\ < ujq. Different quasienergy comes 
near the resonance with the Fermi surface, see Fig. [T] 
The same form of the expression for the tunneling rates, 
Eqs. (|4.5p and (|4.6p . can be used by replacing E Q — > E 
and J m (A) -> J m +,- n {\). 

V. SCATTERING DYNAMICS 

Let us now investigate the modification of tunneling 
rates due to scattering, i.e., the fourth term, V sc , of 
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FIG. 2: Tunneling rates as a function of harmonic field am- 
plitude. The rates are given in terms of zero- field tunneling 
rate amplitude 70. The dotted curve represents the tunneling 
in rate with singular renormalization factor due to scattering 
with a > 0. In this case, only qualitative dependence is pre- 
sented, since the scattering factor will involve regularization, 
as discussed in the text. 

the Hamiltonian (|3.7[) . ignored in the previous section. 
This problem is similar to the problem of x-ray edge 
singularity^ To produce a tractable solution, one has 
to assume a special form of the scattering potential, 
Vkk' —> VukUk' ■ Here, u k is the cutoff function which 
is of the order 0(1) for E k ~ Ep, Xluq and vanishes for 
Ek ~ £00 • We also assume u k to be symmetric. 

One can notice that the scattering is only present when 
the electron leaves the impurity site.— This differenti- 
ates the tunneling in and out processes. To investigate 
the two in a uniform treatment, it is convenient^ to re- 
define the perturbation due to scattering for tunneling 
in as V sc — —Vd^dY kk ,u k u k ic\.c k i. We a ^ the cor " 
responding term to the unperturbed Hamiltonian, i.e., 
Ho — > Ho + J2kk' u kUk'c\ck' ■ This will only redefine 
the noninteracting conduction electron Green's function, 
9kk'{u) = g k {uj)5 kk > - iVY,k" 9k{u)u k u k »g k „ k ,(Lu). For 
tunneling out, one still has V sc = VdS J2 kk , u k u k i c k c k ' . 

To utilize the results of Ref. Ejfj, note that the tunnel- 
ing rates can also be obtained by calculating the time 
derivative of (Td(t)d(t ± 0)S'(oo, —00)). The expression 
correct up to A 2 is 

Hn/out = Re^A 2 J 2 (A) (5.1) 

00 

x J ^'(iTciW^dt^Ocfe^OIX^W^^))- 

—00 

Here, the state |) include the |0) or |1) state of the impu- 
rity for the tunneling in or out cases, respectively. The 
problem is to compute the average 

F(t-t') = Y,u k u k ,(Tcl(t)d(t)d\t')c k ,(t')). (5.2) 

kk' 

Then, the rates are found from F(t) via its Fourier trans- 
form as 

%n/out = Re^ A 2 J 2 (A)F(nw ). (5.3) 

n,k 



With the above redefinition of the scattering perturba- 
tion, the average (|5.2|) can be computed for any mag- 
nitude of the scattering amplitude as a one-body prob- 
lem with time-dependent potential (scattering on im- 
purity), as it was demonstrated by Nozieres and De 
Dominies^ This is possible since we assume that the im- 
purity has no internal degrees of freedom. If one defines 
the times of two tunneling acts by t and t' , the average 
(15. 2|) is found via the time evolution of <p kk i (t, t', t, t') = 
{Tc k (T)c\, (t')) with all the vertices describing the scat- 
tering acts restrained to the interval (t,t'). The overbar 
denotes complete evolution. 

When the region around the Fermi energy is of inter- 
est, \5E\ <C ujq and m = 0, the asymptotic form of the 
scattering can be used^ This adds a singular factor to 
F{oj — > 0), and thus the tunneling rate, of the form 

[Zo/(±6E)} a , (5.4) 

where a — 28/tt — (S/tt) 2 and 5E = Ep — E . In two 
dimensions the phase shift is defined by tan 6 = -kD2V. 
Here, £0 is the cutoff coming from u k . The exponent in 
Eq. (|5.4j) is found within logarithmic accuracy. 

In our case, higher energy terms are present. They do 
not comply with the asymptotic approximation for the 
scattered wave functions used to obtain Eq. (|5.4[) . For 
higher energies, m ^ 0, a short-time dynamics of F(t) 
is necessary. In the case when Vt <C 1, it is possible to 
omit the integral term in Eq. (35) of Ref. [30|. In other 
words, the scattering becomes less important. As the 
result, the corresponding tunneling terms are the same 
as in the previous section up to 0(V/mujo) corrections 
which are negligible provided V/luq <C 1. The tunneling 
due to large V is also clear. The tail of the scattering 
renormalization will be added, as a factor, until approx- 
imately the n ~ V/luq tunneling term. 

Finally, for the case of Eqs. (|4.5p and (|4.6[) and assum- 
ing that \SE\ <C ujq, we obtain 

7in = -2tt A 2 D 2 J 2 (X)9(6E) (5.5) 

and 

lout = - 2ttA 2 D 2 J 2 (X)9(-5E) (j^j (5.6) 

OO 

- 2nA 2 D 2 J2 JlW- 

m— 1 

When the external harmonic field is weak the tunneling 
rates again approach the standard expression (with the 
scattering renormalization). For strong fields, the result 
depends on the scattering exponent as well as on SE, 
unlike in Eqs. (|4.5[) and (|4.6p . The above solution is 
not valid for intermediate values of SE, nevertheless its 
possible form is rather clear and can be inferred from 
Eqs. (|5.5l 15. 6p as far as the external field influence is 
concerned. 
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As it was mentioned in the previous section, the im- 
purity energy level may be well below (or above) the 
2DEG Fermi energy. In this different quasienergy 

enters the resonance with Ep resulting in the singu- 
lar renormalization of the corresponding term. The ob- 
tained result applies with the replacements Eq — > Eq and 
J m (A) -> Jm+m(A), where \E - E F \ < wo- 

The singular factor leads to the two effects, similar to 
those discussed in Ref. [30| for x-ray absorption (emission) 
problem. It destroys the jump in energy dependence for 
tunneling rate if a < 0. When a > 0, the jump becomes 
larger but is finite due to the presence of the spin degrees 
of freedom and temperature broadening of the electron 
density near the Fermi energy which tend to quench the 
singularity. In the a < case, the difference between in 
and out tunneling rates is stronger; the higher tunneling 
in rate (see Fig. [2]) makes "fin/out rates closer near the res- 
onance, except for the values of A, such that Jo (A) = 0, in 
which case the singularity is suppressed. Due to this lat- 
ter fact, one obtains sharp peaks in the low-temperature 
resistivity of 2DEG as a function of A, as will be shown in 
the next section. For both cases of a, the difference be- 
tween "fi n / ou t vanishes for high temperatures, since mul- 
tiphoton transitions to the high energy region become 
possible in both ways. 



VI. CONDUCTIVITY OF 2DEG 



Let us analyze how the strong- field modification of tun- 
neling affects the conductivity of 2DEG. Since the oscil- 
lating electric field is perpendicular to the 2DEG plane, 
it should not influence the conduction electrons directly, 
but only through the interaction with the adjacent im- 
purities. We assume that the impurities are distributed 
with the density Uj, small enough to neglect the inter- 
ference between scattering on different sitesj^L as well as 
the tunneling between the impurity sites. The contribu- 
tion due to other scattering processes are not of interest 
here. They are not affected by the field in our system 
and will add A-independent terms to the total resistivity. 
The conductivity is calculated as a linear response to a 
vanishingly small in-plane dc electric field. 

Conductivity due to impurity scattering is given by 
a — e 2 noT(Ep)/m, where e, no, and m are the elemen- 
tary charge, electron concentration, and effective mass, 
respectively. The scattering time t(Ef) can be esti- 
mated with the Green's function relaxation time. The 
difference between the two is a well known (1 — cos 6*') 
factor.— As far as the field influence is concerned, one 
can estimate r(Ep) by evaluating the imaginary part of 
the retarded self-energy of conduction electrons^ i.e., 

T~ 1 (Ep) — -ImEretC&F, w — * 0). 



For a dilute impurity system^ 

^ S k( iuJ n) = rii \ V k k + - XI V kk'V k 'kGk'(i^n) (6.1) 
I V k' 

+ V ^k'Vk' k"V k "kGk' (Wn)Sk" fan) + •••>, 

V k'k" ) 

where Gkfan) is the Matsubara Green's function of 
2DEG electrons. In the limit of small concentrations, 
rii — > 0, the noninteracting function [without ~E B £(iuj n )] 
may be used. Assuming Vkk' —> Vu^u^ as before, for 2D 
electron system with stationary impurity scatterers, one 
has24lmX^ t (Af>,w-> 0) = -^sin 2 5. 

The tunneling affects the equilibrium occupation of the 
impurity, as well as Gk{iu n ). The scattering vertex is 
modified as Vkk — ► —Vkk{dcft) = — VkkW — P(A)] f° r the 
donor impurity site [for the acceptor site, one has Vkk —> 
— VfefcP(A)]. We note that for the equilibrium state, the 
averages are over |) = V~P(\)\) 1 + \J\ — P(A)|) , where 
|), corresponds to the filled impurity site and |) to the 
empty site. The occupation probability is 

P(A) = ^" (A) j- . (6.2) 

lout(X) +7m(A) 

The 2DEG electron Green's function becomes 
Gkiiuin) = [Go.kii^ny 1 - Efc(iw„)] _1 , where Go,k(iu n ) 
and Y, t k {iuj n ) are noninteracting Green's function of con- 
duction electrons and the corresponding self-energy due 
to the tunneling potential (|4.2p . respectively. The lat- 
ter is Y? k (iu) n ) = J2 m ^ 2j mWG{ioJ n - mu Q ). For the 
impurity scattering self-energy, it is sufficient to use the 
noninteracting function G — > Go m £jUiu; n ). 

When calculating the tunneling contribution to con- 
ductivity, the effect of scattering in S^(iw„) is included 
as suggested earlier. This results in the factor of 
(5E/£o)( s /^ 2 /8E for the term with m = 0, with the 
divergence, at SE — > suppressed, as discussed above. 
The renormalization of A at resonance, i.e., A — > kA, 
takes place. The renormalized term, however, may be 
suppressed by the choice of A, such that Jo (A) = 0. 

Finally, we have two contributions to the conductivity. 
One is due to the tunneling, 

ar 1 ~ rtiEp)- 1 ~ A 2 J 2 (A)<5'(<5£), (6.3) 

where the 5-function 8'(x) is broadened so that i5'(0) = 
k 2 . The other is caused by scattering, 

a- 1 ~ T sc (E F y l (6.4) 

nitan 2 (a) 7 r 2 (A)7 2 Mt (A) 

' 7» 2 (Ah 2 ut (A)(cos<5)- 2 + + f ' 

In Figs. [3] and [H we plot the reciprocal conductivity 
due to tunneling resonance, Eq. (|6.3p . and scattering, 
Eq. (|6.4[) . as a function of A and 1/A. The result is given 
in terms of the conductivity <jq due to the scattering on 
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FIG. 3: Reciprocal conductivity as a function of harmonic 
field strength. The right complementary scale is for accep- 
tor impurity sites. The solid curve is a^c f° r SE > 0. The 
horizontal solid line represents ct^. 1 for SE < 0, when the im- 
purity sites are empty. The dashed curve shows a^ 1 up to 
a factor of ft 2 . In this case, only the left scale is appropri- 
ate. The dash-dotted curve represents ct^. 1 for SE > with 
the singular renormalization factor with a > 0. For stronger 
amplifications, the curve follows closer to zero in its mini- 
mums. The dotted curve gives a qualitative dependence for 
finite temperatures at SE — > 0, a > 0. In this case, the recip- 
rocal conductivity curve becomes bounded from the bottom 
(dotted horizontal line) 




0.2 0.4 0.6 
cj /eE z \z 0i0 - Zk,k\ 



are given by Eqs. (|4.5|1 and (|4.6p . When SE < 0, the im- 
purity sites are empty at equilibrium and the scattering 
occurs with the highest probability — the conductivity is 
not affected by the field. When SE > 0, cr^ 1 oscillates 
as a function of A. For A = 0, the occupied impurities 
do not scatter 2DEG electrons (the situation is opposite 
for acceptor sites). At A 3> 1 tunneling out transitions 
from higher quasienergies dominate, leaving the impurity 
empty. The values of resistivity corresponding to SE > 
and SE < converge to each other, see Fig. @] In the 
figure, tan S = 0.1. For larger values of V the oscillations 
converge to "1" faster. 

When the singular renormalization is introduced, the 
resistivity peaks become sharper for a > and SE > 0. 
The singularity amplifies the tunneling to the impurity 
site for all A except for near the zeros of Jo (A), deactivat- 
ing the scatterers. The corresponding curves in Figs. [3] 
and 2] are given for amplification factor of 10. The esti- 
mate for small (but finite) temperature with SE — )■ is 
shown as a dotted curve in Fig. [3l In this case, the hori- 
zontal line gives the lower bound for large amplifications. 

We have investigated tunneling and conductivity mod- 
ification in impurity-2DEG structure in external time- 
periodic field. Nonlinear dependence on the field ampli- 
tude has been obtained and analyzed for both tunneling 
rates and conductivity. The calculations have been per- 
formed in the limit of small A and nearly zero temper- 
atures. Further investigation has to be done to under- 
stand the modification of other 2DEG electron correla- 
tions, such as the ones leading to RKKY coupling. The 
presence of larger currents is also of interest. 
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FIG. 4: Reciprocal conductivity as a function of harmonic 
field frequency, loq. The right complementary scale is for ac- 
ceptor impurity sites. The curves are cut at low frequency to 
exclude heavy oscillations. The solid curve is a 
The dashed curve is a^ 1 up to a factor of ft 2 
scale is appropriate in this case) ; the dotted curve corresponds 
to the dash-dotted line in Fig. [3] 
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APPENDIX A 



stationary ionized impurities in the absence of external 
harmonic field. The tunneling contribution to resistivity, 
erT" j features the two state aspect of the impurity-2DEG 
coupling. It reflects the dynamical suppression of res- 
onant tunneling (dynamical localization) similar to the 
double quantum dot system o 16 i 17 ' 18 i 27 where the tunnel- 
ing is suppressed by Jo (A) (or Jm(A) for lower Eq or 
biased structures) as well. This contribution is propor- 
tional to A 2 k 2 and vanishes for large fields as 1/A. It is 
independent of the donor (acceptor) type of the impurity. 

The solid curve in Fig. [3] and [4] show the resistivity 
due to scattering of conduction electrons on tunneling- 
active impurities, c" 1 . The modification of tunneling is 
not considered in this case, and the corresponding rates 



We need to analyze the decay of 
(S(-oo,t)d(t)d i (t)S(t,-oo)). The scattering ma- 
trix is due to the interaction (|4.2[) . with the unperturbed 
part (|4.ip . Define the equilibrium Green's function 
G(t,t') = -i{Td{t)tf(t')S(-ao,ao)), Then, the desired 
average is G(t, t' — > t ± 0). 

Since Ho and V are time-independent, one can obtain 
the dynamics by differentiating G(t, t') with respect to 
one of the times. Inverting the operator id/dt' — Eo, one 
obtains 

oo 

G(t - t") = g(t - t") -iAJ2 fdt'M m (t,t')g(t' - t"), 

m ' k -oo 

(A.l) 
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where 

oo 

M m (t,t') = J m (A)^H)" J dh...dt n (A.2) 

x{Td(t)V(t 1 )..y(t n )c{ l (t')I^\t')). 

Here, summation over the connected diagrams is as- 
sumed. It is easily noticed that only odd-n terms will con- 
tribute to the expression. The first-order term contains 
the average (Td(t)I ( l l) (t^ (t 1 )c k (t 1 )cl,(t / )I {m) (t 1 )) 
which splits as g(t — ti)<?fc(ii — t')g m (ti — t'), where 
g(t — ti) and g k (ti — t') are noninteracting Green's 
functions for impurity and 2DEG electrons, respectively, 
while g m (h - t') = e""""^''). The corresponding 
diagram is 



where the sign of summation indices entering the ex- 



where the solid arrow is g(t — t\), the double arrow is 
gk (ti —t'), and the dash-arrow represents g m (t\—t'). The 
next-order, n — 3, terms have two g(t\— t 2 ) andg^ti — t 2 ) 

with four spin vertices I± . Some of the diagrams are 



(a) 



(b) 



The second (b) diagram represents the fact that aver- 
ages (I i + 1 \t 1 )I i I l2 \t 2 )I { + a) (t 3 )I im) (t')) do not necessar- 
ily split in pairs — the only restriction is 72,1—712+713 — m = 
0. The diagrams of type (a) contribute to a single loop 
approximate solution, 



G(uj) = g(uj) + g(w)£i(cj)G(a;), 
where the self-energy is 



(A.3) 



S 1 (w)=^A 2 j2(A) 5fe (w-ma; ). (A.4) 

m,/c 

This result corresponds to the first-order self-energy of 
the general solution. The other terms can also be evalu- 
ated. 

Note that operators I± describe an asymptotically 
large spin2& and thus commute with each other when 
averaged over the states corresponding to zero (finite) 
magnetization, see Eq (|3.9p . Therefore, one can replace 

±77 in the expression equals zero. After some algebra 
the N-ih order (N > 1) self-energy can be obtained in 
the form 

EjvM = Y, A 2N J ni {X)J n2 {-\)... 

7li,...,TO2iV,fci,...,fcjv 

x S(n 1 + n 2 + ... + n 2 N-i + n 2 N) (A. 5) 

x gk 1 ,n 2 .n 3i4 ,... { UJ )9n 3A ,...,n 2 N-l,2N 
X 9n 2N -i.2N{ UJ )9k N ,n 2N {^), 



pression as —77 has been changed, riy = rii 



We 



have defined (cj) = g(ui — loq^ — ujq^' — ...) and 

9k,£,£',...(^) = 9k {u - - ^oC' - ■■■)■ Tne summa- 
tion in Eq. (|A.5[) should not include the terms which 
have g(k)x ...( a; ) with £ + £' + ... = 0. These terms are 
taken into account by the preceding self-energies. An es- 
timate of the above sum suggests that Sat ~ x N , with 
x = A 2 D 2 /u a . In our case, A 2 D 2 /u a ~ A 2 /lu E f <C 1. 
Therefore, a one-loop approximation is sufficient. 



APPENDIX B 

Let us obtain the first-order tunneling rate by averag- 
ing the tunneling amplitude subject to evolution due to 
the Hamiltonian (|2.5[) without the scatterin g te rm. This 
procedure is similar to the one used in Ref. [3fj| to study 
ionization of quantum dot in high-frequency harmonic 
field. The amplitude for an electron to tunnel from im- 
purity site to the 2DEG states is A(d(t)c\.(t)) . Here, 



the overbar denotes evolution due to Hamiltonian (|2.5p 
with only the scattering term absent. The first nonvan- 
ishing order (with the complex prefactor ignored) is 

00 

J dt'Y^(Td(t)cl(t)dHt')d k (t')), (B.l) 



where d(t) = de -iE t-i^nt 



c k e 



"■o . Ihese expressions are substi- 



and c k (t) = 

ie expressi 
tuted into Eq. (|B.1[) to obtain 

00 

f dt'Y A2 { T d(t)ci(t)d^t')c k (t')) (B.2) 

k 

-00 K 

xY J nWJn'We- muJot e m,UJot ', (B.3) 



expanding the oscillatory exponents into the Bessel se- 
ries. Here, d(t) = de- lEot and c k {t) = c k e-' lEkt . The 
amplitude is then averaged over the fast driving field os- 
cillations yielding 



±t nUo t with the restriction that the sum of all J 2 (\) f dt' V A 2 {Td{t)c\(t)d) (t l )c k {t'))e m ^ t ' - f ) . 

expression equals zero. After some algebra, „ J , 



(B-4) 

After the integration, one finally obtains the rates pro- 
portional to 

I m £A^(A)V 6 f ] l n \.^ (B.5) 
E - E k - uj n ± zO 



which gives Eq. (|4.3|) with u) —* Eq 
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